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Abstract. The search for Blank Field Sources (BFS), i.e. X-ray sources without optical counterparts, paves the way to the 
identification of unusual objects in the X-ray sky. Here we present four BFS detected in the Brera Multiscale Wavelet catalog of 
ROSAT HRI observations. This sample has been selected on the basis of source brightness, distance from possible counterparts 
at other wavelengths, point-like shape and good estimate of the X-ray flux (fx). The observed fx and the limiting magnitude of 
the optical catalogs fix a lower limit for our BFS on fx/f op t ~ 40. This value puts them well beyond 90% threshold for usual 
source classes once HRI energy band and proper spectral shape are taken into account, leaving room for speculation on their 
nature. Three BFS show also evidence of a transient behaviour. 

Key words. X-rays: general, Stars: neutron, X-rays: galaxies: clusters, BL Lacertae objects: general, Quasars: general 



1. Introduction 

It is well known that one can discriminate among classes of 
X-ray sources evaluating the ratio between their fluxes in X- 
ray and optical bands (e.g. Maccacaro et al. [T988 , Motch et 
al. Qj998 1 Zickgraf et al. 2003 1, especially if some spectral in- 
formation is supplemented. Blank Field Sources (BFS), i.e. X- 
ray sources with no optical counterparts, stand somehow apart 
from the bulk of X-ray zoology in that they represent rarer, less 
studied objects (e.g. Cagnoni et al. 2002 1. In particular, known 
X-ray sources concentrate at not too high values of fx/fopt, 
e.g., stars have fx/fopt 5=F galaxies and AGN reach fx/fopt of a 
few tens. At larger values of fx/fopt, peculiar populations start 
to show up: some examples are isolated neutron stars, type 2 
quasars, extreme BL Lacs or high-redshift clusters of galaxies. 

The search of BFS aims at objects with high fx/fopt: the 
observed X-ray flux and the upper limit from the absence of an 
optical detection establish a lower limit on the fx/fopt- A careful 
selection of bright BFS is therefore a powerful tool for finding 
population leftovers from usual source classes. The status of 
BFS is somehow a transitory condition that expresses our lack 
of knowledge on the real nature of these sources. New X-ray 
observations can provide temporal and spectral information on 
BFS, while deeper optical/near IR pointings can either reveal 
the counterpart or establish a much stronger lower limit on the 



fx/fopt- With the aid of these information, BFS can be identified 
and classified properly. 



In this paper we present a sample of luminous BFS from the 
ROSAT HRI Brera Multiscale Wavelet catalog (BMW-HRI: 
Panzera et al. 2003 Campana et al. [T9^9l Lazzati et al.[T933}, 
with fx/fopt <: 40 , along with archival X-ray spectral data 
from other observations (when available). This is the first study 
to exploit all ROSAT HRI data for BFS search (see Musso et 
al.[T998 for the first search on a much smaller dataset). A key 
factor in the identification process is the uncertainty of the X- 
ray position (see Cagnoni et al. 2002, Rutledge et al. 2003); on 
this respect, the use of HRI rather than ROSAT PSPC images 
for searching BFS is much superior. 

The paper is organized as follows: in Section 2 we briefly 
describe the BMW-HRI catalog and its cross-identification pro- 
gram with optical, infrared and radio catalogs, in Section 3 we 
introduce the method of sample selection, while the final sam- 
ple of BFS is presented in Section 4. In Section 5 we discuss the 
possible nature of BFS and future search plans are mentioned 
in Section 6. 
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2. The BMW-HRI catalog 

2.1. Catalog description 

The BMW-HRI catalog (Panzera et al. I2UU51 is derived from 
an analysis of the ROSAT HRI data set with a source detection 
algorithm based on the wavelet transform (Lazzati et al. 1999 
Campana et al. fT999l . 

The up-to-date version of the catalog contains 28998 en- 
tries, down to a count rate of ~ 1CT 4 cts s . The total sky 
coverage is 732 deg 2 (~ 1.8% of the sky). For each entry 
name, position, count rate, extension and relative errors are 
provided, along with derived parameters like flux and Galactic 
column density and ancillary information about the pointing. 
Furthermore, results of cross-correlations with GSC2, 2MASS, 
IRAS, and FIRST catalogs are reported (see Section 2.2 and 
TableQ. In the following, some of the basic characteristics of 
the catalog are discussed in order to permit an understanding 
of the BFS selection method. Further details can be found in 
Panzera et al. ( 1277031 . 

Extension. The wavelet transform is particularly power- 
ful in dealing with extended sources. In practice, each de- 
tected source is characterized by its extension, alongside with 
the more usual parameters (count rate and position). Extended 
sources, especially when observed at large off-axis angles, are 
often broken in various point-like sources by other detection 
algorithms. This risk is minimized by the wavelet transform. 

Detection threshold and spurious sources. Each BMW- 
HRI source is characterized also by a detection probability 
value (peak significance). This value is the probability, in units 
of Gaussian cr, that the source is not spurious (i.e. a background 
fluctuation). This probability is the result of numerical simula- 
tions of random fields (without sources), repeated for different 
wavelet transform scales; it takes into account the number of 
background peaks mismatched as sources by the detection al- 
gorithm, for a given signal to noise ratio threshold in wavelet 
space. The mean number of spurious sources expected over an 
image is kept constant (0.4 spurious detections per field), vary- 
ing the detection threshold. That is, the detection threshold de- 
pends on the mean background counts value on the image and 
on the scale at which the search is performed (see Lazzati et 
al. QP99). The mean threshold value in the overall catalog is 
~ 4.2cr. 

Errors. In principle, the error for each parameter found 
with the wavelet transform algorithm can be estimated from 
the covariance matrix. In practice, such errors are reliable 
only when the number of source and background counts are 
sufficiently high (> 2 x lCT" counts per pixel). Otherwise, 
the distribution of the wavelet transform coefficients becomes 
Poissonian rather than Gaussian (in the lowest scales). If this is 
the case, a better estimate of the errors can be given by means 
of basic statistics (such as the error on the number of counts 
N is yN, for more details see Lazzati et al. I1999i . The errors 
reported here and in the catalog are the maximum between the 
covariance matrix and the statistical estimates by Lazzati et al. 
(19991. 

Boresight correction uncertainty. There is an additional 
source of errors on the absolute position determination of HRI 



sources, the uncertainty on the boresight correction. If the 
alignment between the node and the telescope optical axis is 
perfect, on-axis images are exactly in the center of the instru- 
ment. However, in the real situation there is always a finite mis- 
alignment that needs to be corrected. Unfortunately, uncertain- 
ties in the aspect solution (that describes the orientation of the 
telescope as a function of time) and errors in the alignment be- 
tween the star trackers and the ROSAT X-ray Telescope intro- 
duce uncertainties in the boresight correction, of variable size 
for each observation. In practice, the systematic offset between 
accurately known optical positions and X-ray positions can be 
used to evaluate the extent of boresight correction uncertain- 
ties. For the ROSAT HRI, the offset can be as large as 10", 
even if in most cases it will be much less (David et al. 1998). 
The 10" value is conservatively assumed as a fiducial value for 
the boresight offset. Usually, the statistical errors on the posi- 
tion from the detection algorithm, calculated as in the previous 
section, are much less than 10" and therefore can be neglected 
in a first approximation. 

Multiple detections. As each pointing is treated separately, 
not all of the catalog entries correspond to independent sources. 
An estimate of the overall number of independent sources can 
be given compressing the catalog in a 10" radius (again, the 
fiducial error value indicated by the boresight uncertainty). 
This procedure selects only a source for each 10" cone ra- 
dius, on the basis of an autocorrelation of the position. With 
this procedure, 20433 sources are left. Obviously, the compres- 
sion brings to the loss of sources truly close to each other. 

2.2. Multiwavelength catalog cross correlations 

In the following, we will largely use the term off-band to in- 
dicate a wavelength passband different from X-rays. The lim- 
iting flux of a given catalog in a given band will be generi- 
cally denominated foffband- For each of the BMW-HRI entries, 
cross correlations with GSC2 (McLean et al. I2005t . 2MASS 
Second Data Release (Kleinmann et al. IT9"9"4T l. IRAS Point 
Source Catalog (Beichman et al. [T9881 and FIRST Survey 
Catalog (White et al. [1997) were performed by Panzera et al. 
(2003 ). Off-band catalog properties are summarized in Tabled 
The adopted radius for the cross-correlation is 10" (see Section 
2.1), assumed as positional X-ray uncertainty. In fact, the posi- 
tional uncertainties for GSC2, 2MASS SDR and FIRST cata- 
logs are < 0.5" (3<x, McLean et al. l2005> . < 0.5" (lcr, see the 
2MASS documentation for an extensive discussion), < 0.5" 
(90%, McMahon et al. 2001 1 respectively, so that they can be 
safely neglected for cross-identification purposes. An excep- 
tion is the case of IRAS PSC, for which the positional accuracy 
varies with source size, brightness and spectral shape and it is 
different in different directions, but it is usually better than 20" 
(see Beichman et al. [T9~88 1, so this last value has been used as 
cross-correlation radius with this catalog. Note that: a) when 
two or more entries in the correlating catalog are found, only 
the brightest source parameters are reported (but in any case 
the number of found sources is given), b) the sky coverage of 
the off-band catalogs is usually not complete (see again Table 
[0. In particular, the preliminary, unpublished version of the 
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GSC2 (GSC2.3) catalog used in Panzera et al. j2"UU5t lacks 
coverage of the zones where a bright source caused an over- 
exposure of the Schmidt plates. For the 2MASS and FIRST, 
more complete catalogs are now available (2MASS All Sky 
Data Release, FIRST Survey Catalog 03Aprl 1 Version). Our 
final list has been checked with them. 



Table 1. BMW-HRI cross-identifications. 



Catalog 


Survey 


Band 


Coverage 


Depth 


GSC2 


POSSII 


BJ 


Dec: > 


22.5 


(prel.) 




R 


Dec: > 


20.8 






I 


Dec: > 0* 


19.5 




SERC 


BJ 


Dec: < 


23 




" +AAO 


R 


Dec: < 


22 






I 


Dec: < 0* 


19.5 


2MASS SDR 


2MASS 


J 


50% sky 


15.8 


(2000) 




H 


50% sky 


15.1 






K 


50% sky 


14.3 


IRAS PSC 


IRAS 


12 /im 


98% sky 


0.4 Jy 


(1989) 




25 fim 


98% sky 


0.5 Jy 






60 ytim 


98% sky 


0.6 Jy 






100 yum 


98% sky 


1.0 Jy 


FIRST SC 


FIRST 


1.4 GHz 


20% sky 


~ 1 mJy 


(2001) 











* work in progress. 



3. Sample selection 

3.1. Catalog mask 

The first sample selection was made on the basis of the follow- 
ing criteria: 

a. blank field sources, i.e. without any cross identification in 
the other databases (see Section 2.2 and Table^: 6061 cat- 
alog entries; 

b. point-like sources (significance of the extension=0): 4955 
catalog entries; 

c. bright sources, i.e. with fx > 2.7 X 10~ 13 erg s _1 cm~ 2 : 275 
catalog entries. 

The reported numbers are catalog entries, i.e. no compression 
was applied to eliminate multiple detections (see Section 3.2 
below). The fx we used is computed from the observed source 
counts, in the 0.5-2 keV range, considering only channels 2-9 
(e.g. David et al. 1998 1. The assumed spectral shape is a power 
law with photon index F=2 (i.e. the same of the Crab pulsar 
and nebula) and column density 5 x 10 iy atoms cirT 2 (param- 
eter # 48 in the BMW-HRI catalog). The chosen model for 
the count rate-flux conversion uses constant negligible column 
density, instead of the integrated Galactic one in the source di- 
rection. While the latter is commonly adopted in extragalactic 
surveys, it would introduce a larger error on low-latitude galac- 
tic BFS, which are expected to be low-luminosity, nearby ob- 
jects (see Section 5). Furthermore, the chosen fx derived with 
this choice is lower than the one with Galactic Nh, so the ob- 
tained fx/foffband limit is lower. We chose to select sources on 



the basis of fx and not of, for example, fx/fBj- This brings nat- 
urally to different fx/foffband ratios for sources in different parts 
of the sky, due to the variable depths of the surveys available to 
us. However, in such a way the selection is based on a directly 
observed quantity, the number of counts, not on the basis of 
a non-detection, and it provides a more homogeneous sample 
for future observations. We carefully evaluated the fx/foffband 
keeping into account the different depths and photometric band 
used in the different catalogs. The obtained lower limit for 
the fx/fopt from the GSC2 catalog are fx/fBj > 37 (Northern 
sky), fx/fBj > 59 (Southern sky), f x /fRF > 43 (Northern sky), 
fx/fRF > 129 (Southern sky). A detailed description of method 
and used values can be found in Chieregato (2005 1. 

In conclusion, at odds with the majority of similar stud- 
ies, our sources are blank at various wavelengths (see again 
Table [Q. We note also that, although we selected point-like 
sources, we cannot exclude that the real nature of some source 
is extended (i.e. only the brightest peak has been detected with 
HRI). 

3.2. Multiple detections removal 

The removal of multiple detections from the sample could be 
done, in principle, using automatic compression criteria (posi- 
tion based), like in Section 2.1. However, given the limited en- 
tity of our sample of BFS candidates (275 elements), we pre- 
ferred to remove multiple detections by hand to avoid, when 
possible, the loss of sources truly close to each other. 

We used positional coincidence as the main criterion for 
identifying multiple detections. Another useful piece of in- 
formation was the flux level. Unfortunately, positional coin- 
cidence can be weakened by the boresight uncertainty, while 
obviously the flux criterion is not fulfilled for highly variable 
sources. In practice multiple detections were reasonably iden- 
tified in most of the cases. After multiple detections removal, 
we are left with 226 catalog entries. 

3.3. Total counts cut 

We made a further screening on the number of total counts, re- 
jecting sources with less than 25 photons (on the original detec- 
tion, i.e. before applying vignetting and Point Spread Function 
loss corrections). This has been done in order to strengthen the 
derived source parameters (obviously, due to the statistics in- 
crease). After the total counts cut, we are left with 201 sources. 

3.4. Are the BFS sources real? 

We re-evaluated the number of spurious sources expected in 
this 201 subsample, again on the basis of simulations of ran- 
dom fields with the above masks. In particular, the flux cut is 
very effective in the removal of spurious sources. In fact, the 
expected number of spurious source is ~ 2 in the 201 subsam- 
ple (see Moretti et al. 2004 1. Our BFS are therefore likely not 
to be fake detections. 
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3.5. Visual inspection: when a blank is not a blank 

We inspected Digitized Sky Survey 2 (DSS2; McLean et al. 
20001, Super Cosmos Sky (SSS; Hambly et al. EOOTl . 2MASS, 
FIRST, Sy dney U niversity Molonglo Sky Survey (SUMSS; 
Bock et al. [T9^9t . NRAO VLA Sky Survey (NVSS; Condon 
et al. [T998l an d Wes terbork Norther Sky Survey (WENSS; 
Rengelink et al.[T997i images for each of the 201 sources left. 
Furthermore, we cross-checked source positions with Simbad, 
NED and with the VIZIER catalogs (in particular the USNO- 
Bl catalog, Monet et al. 20031. We excluded from the sam- 
ple sources and fields already known and well studied or too 
complex (e.g. Magellanic Clouds, M31, Orion Nebula). We 
excluded also sources positionally coincident with bright stars 
(holes in the GSC2 coverage), or at the periphery of optically 
extended emission (more than 10" from the centroid and so 
escaped the automatic cross-identification). We found, and ex- 
cluded, some cases in which our X-ray sources coincided with 
reliable USNOB1.0 entries (the criteria for the GSC2 detection 
are somewhat more conservative), or with counterparts at other 
wavelengths (IRAS Faint Sources Catalog, Moshiret al. 1990 1. 
We prefer not to be too severe at this stage, deferring strict cuts 
after the boresight correction. 72 sources survived this phase. 

3.6. Cheshire cat: elimination of fluctuations 

We put aside from the primary sample sources which have not 
been detected in a longer HRI or ROSAT PSPC pointing (the 
PSPC had the same passband but a substantially larger effec- 
tive area than the HRI). These could be variable sources de- 
tected at peak flux, but they are also candidates for false detec- 
tions. No source has been inserted in this subclass on a non- 
detection basis with other instruments, because of the different 
passband, i.e., given the absence of spectral information on our 
sources, the non-detection could be ought to a peculiar spec- 
tral shape. For some sources the suspicion of a false detection 
is strengthened by the number of non-detections or by the re- 
spective length of the pointings. However, there are cases for 
which source parameters like detection probability and number 
of counts seem to indicate a real source. Furthermore, as dis- 
cussed in Section 3.4, we expected only ~ 2 spurious detection 
in 201 sources. Of the 72 sources, 16 were placed in this tran- 
sient candidate sub-sample, while for the other 56 there was no 
evidence of a transient nature. 

3. 7. Pin-pointing the sources to the sky: boresight 
correction 

We performed some astrometry on the remaining fields (includ- 
ing those of transient candidates), in order to get rid of the 10" 
fiducial boresight uncertainty and therefore to fully exploit the 
HRI angular resolution capability. Furthermore, the error given 
by the detection algorithm (see Section 2.1), while not being 
entirely statistical, can be as a first approximation treated as 
random, so that the usual Gaussian relations can be used. In 
contrast, the offset given by the boresight uncertainty varies 
in a random way in the ensemble of all the pointings, but is 



systematic in nature for all the sources in the single pointing, 
weakening the rejection of distant optical associations to BFS. 

For each pointing, we matched X-ray sources to optical 
counterparts (optical positions of known X-ray emitters), if 
any, or to optical catalog sources distant less than 10". We 
excluded X-ray extended sources and sources with no optical 
catalog entry in 10". Ambiguous cases, i.e. with two or more 
optical sources present, were treated individually using the dis- 
tance and the optical luminosity as criteria for the identifica- 
tion and then checked a posteriori. Even if only another X-ray 
source was present, the shift for this source has been applied 
to the pointing, although obviously in these cases the corrected 
positions and uncertainties have to be taken cum grano salis. 
Generally, the new error is less than 10"; even if it remains 
around 10", the new position should be free from systematic 
boresight uncertainty. Note that performing the boresight cor- 
rection only after the selection of Section 3.1 can bring to the 
loss of BFS. However, this allows us to deal with a limited 
number of sources and boresight corrections. 

4. The data 

The elimination of sources with an off-band counterpart in a 
4<x radius from the boresight corrected position left us with 
our final sample of four sources (three of them transient candi- 
dates). This Act limit assures us that starting with 1340 sources 
(i.e. sources outside the inner 3' in the BMW-HRI catalog and 
obeying to the cuts above unless th lack of counterparts), only 
0.08 of them would lie outside the error region by chance. 
In summary, our final BMW-HRI Blank Fields Sources are: 

- Blank: no counterpart at other wavelengths in a 4<x radius 
from the boresight corrected position. 

- Bright: f x > 2.7 x 10~ 13 erg s _1 cirT 2 . Note that from the 
GSC2 magnitude limit, this is equivalent to fx/fBj respec- 
tively > 37 and > 59 for Northern and Southern sources. 

- Point-like: extension significance is 0. 

- Well detected: total (uncorrected) counts > 25. 

After this selection we end up with one persistent BFS plus 
three transient BFS. We report the new positions and error radii 
(as obtained from the boresight correction procedure), the flux, 
the detection probability (i.e. the probability to detect a back- 
ground fluctuation as a source), the total number of counts, the 
Galactic coordinates, the integrated Galactic column density, 
the distance from the nearest off-band association (in terms of 
error radii), the lower limit on the fx/fBj and the upper limit on 
the radio emission (see Table|2}. Note that, though with differ- 
ent depth, all but one (lBMW200739.8-484819)of the sources 
have a radio flux upper limit. 

4.1. Properties of BFS 

The final error radius estimate for our sources is quite different 
from case to case, ranging from ~ 3" to ~ 7"; the latter large 
error is due to the large off-axis angle and to the low number of 
counts (see Fig.[Q. In fact, all our sources are quite offset from 
the center of the pointing (> 15')- This depends from the used 



Table 2. BFS parameters. 
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Name 


R.A. 


Dec 


bll 


Err. rad. 


Flux 


Prob. 


Cts 


N fl 


fx/f B J 


Close ass. 




J2000 


J2000 




tf 


erg cirT 2 s _I 


cr 




10 20 cm- 2 




(sigma) 


1BMW042142.4-571541 


04 21 42.6 


-57 15 39 


-42.54 


2.9 


6.5x10-" 


14.0 


742 


1.8 


> 141 


5.3 


1BMW135703.0+181122 


13 57 02.9 


+ 18 11 21 


72.45 


7.1 


3.5x10-" 


4.2 


112 


2.1 


> 47 


6.3 


1BMW200739.8-484819 


20 07 39.8 


-48 48 19 


-32.31 


5.8 


3.0x10-" 


4.3 


55 


5.1 


> 65 


5.2 


1BMW043306.8+ 155307 


04 33 06.8 


+ 15 53 07 


-21.16 


6.5 


2.9x10-" 


4.2 


34 


17.0 


> 40 


4.9 






5 10 15 20 

off axis angle [arcmin] 

Fig. 1. BMW-HRI position error: the statistical errors on posi- 
tion of BMW-HRI sources for different number of counts as a 
function of the off-axis angle. Error bars are at 2<x. The change 
in the slope of the curves corresponds to angular resolution loss 
at the change of image rebinning. An offset of 0.2' and 0.4' has 
been applied to points in the 50 counts bin and 100 counts bin, 
respectively, for visualization purposes. 

selection procedure and needs some comments. First, we con- 
centrate on the consequences of the requirement of absence of 
offband counterparts in a 10" radius. This affects strongly the 
overall distribution in detector coordinates. In fact, the fidu- 
cial 10" radius for the initial correlation, chosen to match the 
worst case boresight error, is huge in comparison to the sta- 
tistical positional error for the innermost sources, but it cor- 
responds to the intrinsic error for the outermost sources (see 
Fig-0. Sources in the innermost part are preferentially associ- 
ated with optical counterparts which may not be the true ones, 
since they have been searched for over a region much larger 
than their error boxes. The overall effect is to deplete the dis- 
tribution from inner sources, enhancing the proportion of outer 
sources. This bias is present independently of the real size of 
the boresight correction error. Second, there is a trend in over- 
estimating outer source counts (especially at low values), since 
it is favorable to detect them on top of a positive background 



fluctuation (an effect known as Eddington bias, see Hasinger 
et al. 119931 Moretti et al. 2002 1. This again adds on the en- 
hancement of outer sources. The first bias leads to the loss of 
BFS, but it does not affect our selected sources. The effect of 
the second bias on the estimated flux can be evaluated a pos- 
teriori from the duration of the pointings and the number of 
counts. Very conservative assumptions give a factor of ~ 2 as 
the maximum overestimate. 

SSS btat 



1BMWW21A2 4-571 6*1 



o 




10" 



Fig. 2. The X-ray position of 1BMW042142.4-571541 super- 
imposed to the corresponding blue image along with it lcr error 
circles. 



4.2. 1BMW0421 42.4-571 541 

This source is the brightest of our sample and the one with 
the highest number of detected photons. The good statistics in- 
volves a small detection algorithm error radius, so the position 
determination is quite accurate (see Hig. EJ. The boresight cor- 
rection was performed with the known position of the quasar 
HE 0419-5657 and the nucleus of NGC 1547 as well as with 
other eleven sources matched to optical catalogs. The source 
is close (~ 15" but actually more than 5 cr) to a bright star 
(BJ=13.31, F=11.61 in the GSC2.2); this makes less reliable 
the lower limit on the fx/f pti m particular in the red band. 
There are two fake USNO B1.0 sources along the saturation 
spikes of the bright star (see Fig.|3}. The source was observed 
with Einstein (IE 0420.7-5723) and with the ROSAT PSPC 
(1WGAJ042 1.7-57 16) and was also detected as a RASS Bright 
Source (1RXS J042144.0-571601). The WGA computed flux 
is 8.8 x 1Q-' 3 erg cm -2 s _1 , consistent with the RASS count 
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rate. Taking into account the caveats on the count rate to flux 
conversion, this could mean that there is no evidence for long 
term flux variations. 

We re-extracted the archival PSPC observation of this 
source (sequence rp700034n00). The source was very close 
to the PSPC rib; despite this, we attempted a rough spectral 
analysis. We extracted source photons from a circle of radius 
200" centered on source position. We used as background re- 
gion a circle of radius 240", centered in an empty region south 
of the source. The collected counts were 321. We rebinned 
channels by a variable factor in order to have at least 20 pho- 
tons in each bin. Channels 1-11 and 136-256 were ignored. 
Figure 3 shows the resulting spectrum, which we then fitted 
with XSPEC (v. 1 1.2). Due to the poor statistics, we considered 
a limited number of models, as shown in Table [3] The spec- 
trum is definitely soft, without evidence of absorption excess. 
The best fit absorbed blackbody model tends to null absorption, 
with a 90% upper limit on column density of ~ 3 x 10 19 atoms 
cm" 2 . The power law fit gives a relatively better reduced x 2 
(1.3 versus 1.4). Furthermore, the best column density value is 
quite close to the integrated Galactic value of 1.2 X 10 20 atoms 
cm" 2 (see Fig.[3J- Freezing absorption at the Galactic value, the 
reduced x 1 does not change appreciably. The bremsstrahlung 
fit gives the best reduced^ 2 (1.3) and the column density is 
~ 1.3 x 10 20 atoms cm" 2 . 

Table 3. Spectral fits for PSPC observation of 
1BMW042142.4-571541. 



Model 


N fl 


kT /T 


Flux(0.1-2keV) 


^red 




10 2 " cm" 2 


keV 


10" " cgs 




B.body 




19 +U02 


5.6 


1.4 


B.body 


1.8 (frozen) 


15 +002 
u - 1J -o.oi 


6.9 


1.9 


Pow.-law 


lJCJi 




6.5 


1.3 


Pow.-law 


1.8(frozen) 


1.8- 2 


6.5 


1.2 


Bremss. 


1 3 +1 ° 


1 46 +14 


6.2 


1.3 


Bremms. 


1.8 (frozen) 


1 

LAI3 -0.25 


6.2 


1.2 



4.3. Transient BFS 

In this section are presented three sources (see Table [2} that, 
while fulfilling all the others criteria for BFS, show additional 
strong evidence of a transient behaviour. In fact, despite their 
brightness in the BMW-HRI catalog (f x > 2.7 x 10" 13 erg s" 1 
cm" 2 ), for each of them there is at least another ROSAT HRI or 
PSPC pointing where they have not been detected (deeper than 
the HRI exposure that revealed them). 

1BMW135703.0+181121. 1BMW135703.0+181121 has 
been detected in a ~ 5.5 ks HRI pointing in Jan. 1998; we in- 
cluded it in the possible transient sub-sample since it has not 
been detected in three longer PSPC pointings (~1 1 .4 ks, Jan. 
1992; ~8 ks, Jul. 1990; -5.6 ks, Jul. 1992; also it was not seen 
in a shorter ~1.1 ks HRI pointing in Jan. 1992). The high num- 
ber of counts detected points toward a real, transient source. 
The boresight corrected position has been obtained with three 
sources. There are no counterpart at other wavelength in a 45" 
radius (~ 6.3cr, see Fig.0J. 
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Fig. 4. The X-ray position of 1BMW042142.4-571541 super- 
imposed to the corresponding blue image along with it lcr error 
circles. 
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Fig. 3. Power law spectral fit of the PSPC observation of 
1BMW0421.4-571541. 



1BMW200739.8-484819. This source has been detected in 
a ~ 3.2 ks pointing in Oct. 1996, while no corresponding source 
was observed in four other pointings of comparable length (2.6- 
3.6 ks) at distance of ~ days, nor in two ~1 1 ks ROSAT PSPC 
pointings in Nov. and Apr. 1992. The closest optical catalog 
source is at more than 30" (see Fig.|5j from the boresight cor- 
rected X-ray position (the boresight correction has been per- 
formed with the other two X-ray sources detected in the Oct. 
1996 HRI pointing). Two objects, A and B in Fig.|5] are clearly 
visible in the blue SSS image, respectively at ~ 16" (~ 2.7<x) 
and ~ 13" (~ 2.3cr). Object B (BJ=22) is reported to be ex- 
tended, object A (BJ=22.3) is possibly extended too. Oddly, 
while nothing is apparent at that position in the red SSS im- 
age, a source can be seen both in the shallower ESOI RED 
image and in the SSS I image, at a position almost coincident 
with object B, with RF=20.2 (more than 2 mag brighter than 
in the Schmidt plate limit) and IN=19.0. An obvious hypoth- 
esis could be the presence of an optically variable source, that 
perhaps could be identified with 1BMW200739.8-484819. 
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Fig. 5. The X-ray position 1BMW200739.8-484819 superim- 
posed to the corresponding blue image along with it 1 <x error 
circles. 

1BMW043306.8+155307. The observation in which this 
source was detected (~2 ks in Feb. 1998) is part of the Hyades 
ROSAT campaign. The source is included among transients 
due to its non-detection in four PSPC pointings (of ~2A ks, 
Feb. 1991; -1.9 ks, Mar. 1991; -16.6 ks, Sep. 1992; -19.6 ks, 
Aug. 1993). In all the four PSPC pointings the source was at 
much larger off-axis angle than in the HRI one (31.5' - 55.1' 
versus 15.8'). The boresight correction was done with three 
sources. There are no possible counterparts in a ~ 25" (~ 3.8cr) 
radius from the source position, where there is a source de- 
tected only in the APM North Catalog (E=19.8), not visible in 
the DSS2 digitized sky images. The other closest source is at 
32"(~4.9cr, see Fig. iJ. 




Fig. 6. The X-ray position 1BMW043306.8+155307 superim- 
posed to the corresponding blue image along with it lcr error 
circles. 



5. Discussion: nature of the blank fields 

In order to proceed further with hypotheses on BFS nature we 
need to compare the fx/fBj limits of BFS with existing classi- 



fication schemes. We chose as a reference scheme the one of 
the RASS-Hamburg optical identification program (Zickgraf et 
al. 2003), which uses similar X-ray (0.1-2.4 keV) and optical 
(Johnson B) bands. We converted the RASS-Hamburg limits 
(see Table 6 of Zickgraf et al. l2003l to BMW-HRI and GSC2 
fx/fBj (the deepest GSC2 band), keeping into account the dif- 
ferent X-ray energy band and spectral shape used for flux eval- 
uation, and the different optical band. The resulting classifica- 
tion scheme is shown in Table |4] (see Chieregato 2005 for de- 
tails about the conversion). The fx/fBj lower limits of BFS, ^ 
37 in the Northern sky and ^ 59 and in the Southern sky (due to 
different GSC2 depth, see Tabled, are a factor of ~ 2 and ~ 3 
beyond the limits of our reference values for different classes, 
respectively. This points out the possibility that BFS are un- 
usual objects, left outside from classifications, or at least very 
peculiar members of the more ordinary categories. 

Table 4. fx/ fsj for various classes of sources, adapted to 
BMW-HRI fluxes and GSC2 BJ band. 



Class 


Lower limit 


Upper limit 


B-V 


Stars 


1.3xl0- 5 


0.6 


-0.5-2 


White Dwarfs 


1.2X10" 5 


9.3xl0~ 4 


-0.3-1.5 


Cat. var. 


i.8xitr 2 


4.6 


-0.1-1 


Galaxies 


3.ixi<r 2 


11.1 


0.4-1.5 


Gal. clusters 


0.2 


16.9 


1.0-1.5 


AGN (w. BLLacs) 


0.22 


12.0 


-0.5-1 



5.1. BFS as unusual sources 

It is well known that some rare class of sources can reach ex- 
treme values of fx/f pt (see Cagnoni et al. 2002 1. Here we will 
focus on our specific sample and discuss the possibility that our 
BFS belong to the following categories: 

a. Isolated Neutron Stars (INS). 

b. X-ray binaries. 

c. High Redshift or Dark Clusters of Galaxies. 

d. Type 2 Quasars. 

e. Extreme BL Lac. 

Isolated Neutron Stars. Isolated Neutron Stars (INS; see 
Treves et al. 2OO0| Haberl EUOll Haberl l2DM are extreme BFS. 
The bona-fide fx/f pt lower limit used to assess their nature is 
1000, but for optically identified objects the real value can be as 
high as 10 5 . The X-ray spectrum is optimally fitted by a black 
body with kT~60-100 eV. The optical emission lies a factor 
of a few over the Rayleigh- Jeans tail of the blackbody and no 
emission outside the optical/soft X-ray range has been found. 
The low column densities from the X-ray spectra and the paral- 
lax measured distance for RXJ 1856.5-3754 (Walter & Lattimer 
2002 1 makes this source an intrinsecally faint and closeby ob- 
ject, probably emitting from the neutron star surface. 
Different kind of isolated neutron stars can also reach 
high fx/fopt values: Anomalous X-ray Pulsars (AXP; Israel, 
Mereghetti & Stella 2002 1, close-by radio pulsars with sub- 
stantial cooling X-ray thermal emission (e.g. PSR 0656+14) 
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and Geminga-like objects (Bignami & Caraveo ll996l l. AXPs 
are young, luminous and distant objects and should be much 
closer to the Galactic plane than our BFS. The radio limit, 
though not strong, argue against close-by cooling radio pulsars. 
This does not apply to Geminga-like neutron stars for which 
the radio emission is only marginal. Converting the X-ray flux 
of Geminga to the HRI pass-band and using V=25.6, yields 
fx/fv ~ 67, making it concistent with our BFS candidates. In 
the cooling framework, the same physical mechanism is at the 
basis of INS emission and of the soft thermal component of 
young radio pulsars and Geminga-like objects. Basing on this 
scenario, Popov et al. (2003, 2005 1 produced the log N-log S 
of cooling neutron stars (all the three categories), allowing us 
to estimate that only ~ 1 unidentified INS is expected for our 
flux and sky coverage. If our steady BFS (1BMW042142.4- 
57 1541)) is an INS, no optical counterpart could be revealed up 
to the ~ 28 mag; therefore, this possibility could be strength- 
ened by the absence of optical sources in deeper observations, 
reaching the critical fx/fopt lower limit of 1000, correspond- 
ing to magnitude ~ 25-26, together with a better position from 
Chandra. At the same time the well-defined properties of the 
X-ray emission could be tested by observations with an in- 
strument with sufficient soft response, like the XMM-Newton 
EPIC pn camera. 

X-ray binaries. Some kind of X-ray binaries can reach 
high values of fx/f pt- In particular, BFS searches have already 
discovered Ultraluminous X-ray Sources (ULX; e.g. Cagnoni 
et al. 2003). However, the observed flux level of BFS is too 
low for a Galactic object and none of our BFS can belong to 
nearby galaxies, excluding ULXs, supersoft sources and burst- 
ing LMXBs. Therefore, we consider unlikely known classes of 
X-ray binaries (see however the discussion on possible tran- 
sients). 

Clusters of Galaxies. BFS research has already been ef- 
fective for the discovery of clusters of galaxies since Cagnoni 
et al. (2001 ). Furthermore, several distant clusters have already 
been discovered in the BMW-HRI catalog as extended objects 
(Moretti et al. 2004). Distant clusters may escape this selection 
if only the X-ray peak is bright enough to be detected, resulting 
in a point source, while the rest is concealed in the background. 
This effect is strongly favoured at large off-axis angle (and con- 
sequent PSF degradation) like those at which BFS are detected. 
Therefore, if some BFSs are galaxy clusters, we expect them to 
be at quite high redshift. The fx/f pt of clusters of galaxies is 
the highest among the X-ray emitting classes of Table |4] For 
normal galaxy clusters, the values reached are usually lower 
than those of BFS (and if only the peak of the X-ray emission 
is computed, the resulting fx/f pt is further lowered). Despite 
this, there are two effects that could boost the fx/f pt- The first, 
and more effective, is the redshift of the optical emission. In 
fact, for a redshift of 0.75, the 4000 A break of the cD galaxy 
(or of any early type galaxy of the cluster) would be redder than 
the GSC2 F band. Even a lower redshift could result in a suffi- 
cient enhancement of the fx/f pt- Second, the X-ray peak could 
be not coincident with the optical one (if the cluster is not viri- 
alized, or if there is no cD galaxy, like in Bautz-Morgan type 
III clusters, or if there is a central cooling flow). Optical-near 
IR deeper observations should reveal a substantial increase of 



galaxy counts in proximity of BFS and, if there is any, the cD 
galaxy of the cluster. An even more strong test would be X- 
ray observations with Chandra or XMM-Newton, that should 
reveal an extended source. 

BL Lacs. BL Lacs stand somehow apart from the bulk of 
AGNs in that they can reach high fx/f pt values. The unified 
spectral energy distribution of BL Lacs (Fossati et al. [1998, 
Ghisellini et al. 1998 ) can be described by a two peak emission, 
in which the lower energy peak (radio to X-ray) is commonly 
ascribed to synchrotron emission, while the high energy peak 
(X-ray to y-ray) is probably due to inverse Compton emission. 
In this scenario, different BL Lacs are distinguishables on the 
basis of the shifts of the peaks, creating the so-called blazar 
sequence. By tuning the position of the peaks, it is possible to 
produce the old categories of X-ray selected and radio selected 
BL Lacs, the former being those capable of high fx/f pt values. 

In order to check the viability of BL Lacs hypothesis for 
BFS we consider the different a xo (the power law index con- 
necting the X-ray flux to the optical flux) and a xr (the power 
law index connecting the X-ray flux to the radio flux) for 
Einstein Slew Survey Bl Lacs (Perlman et al. 1996) and the 
upper limit for BFS. Despite the poor radio limits, the locus 
occupied by BFS is partially superimposed to the most ex- 
treme Slew BL Lacs, with the brightest BFS, 1BMW042142.4- 
571541, exhibiting an a xo limit substantially lower. If these ob- 
jects were indeed BL Lacs, they would probably stretch the 
sequence in the a xo -a xr plane, representing a still undiscovered 
extreme population. The identification of the optical and possi- 
bly of the radio counterpart would be very important for testing 
the BL Lac hypothesis. In this case, further X-ray observations 
should show a power law spectrum with photon, without evi- 
dence of absorption excess. 

Type 2 Quasars. Evidences for a substantial population 
of type 2 Quasars, as requested by unification models (Urry 
& Padovani 1995; Comastri et al. 1995), has been continu- 
ously increasing in the recent past. In particular, the resolu- 
tion of the largest fraction of the cosmic X-ray background in 
discrete sources (Hasingeret al. 1998 Mushotzky et al. 2000, 
Campana et al. l20"0"Tl Giacconi et al. 120021 Moretti et al. 120031 
Brandt & Hasinger 2005 ) has brought to infer the existence of 
an adequate number of highly absorbed and luminous AGNs, 
responsible for the hard part of the background. Up to now, 
while only a few type 2 quasars have been firmly identified 
(e.g. Stern et al. 2002; Norman et al. 2003), many candidates 
are being produced either by deep fields, pencil beam searches 
(Chandra Deep Field North, Barger et al. 120031 Chandra Deep 
Field South, Szokoly et al. 2003 Lockman Hole, Mainieri et 
al. [2002), or by dedicated shallower surveys, with a substan- 
tial sky coverage (Hellas2XMM, Baldi et al. I^UUl ChaMP, 
Green et al l2003l SEX SI, Harrison et al. r2003l BMW-Chandra, 
Romano et al. 2005 ). The most striking characteristic of type 
2 quasars is the severe absorption excess. In particular, this 
makes them Extremely Red Objects, with R-K> 5. Therefore, 
if some BFS are type 2 quasar, we expect them to have a rel- 
atively bright near-IR counterpart. Unfortunately, the near-IR 
limit of the 2MASS (Tabled is not sufficient to confirm or 
rule out this possibility. Deeper near-IR observations are there- 
fore mandatory to test the viability of this scenario. Future X- 
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ray spectral information could also be extremely useful, con- 
firming or excluding the need of absorption excess to fit the 
spectra. A number of source capable of very high fx/f pt ratios 
was recently discovered in Chandra Deep Fields, the so-called 
Extreme X-ray Objects (EXO; Koekemoer et al. 2004 ). The na- 
ture of these sources is still not clear, but they share with type 2 
quasars the extremely red colour. However, these sources have 
X-ray fluxes much fainter than BFS (< 4 x 10~ 15 erg cirT 2 s~'), 
so they are likely to belong to different populations. 

The nature of BFS reported here is probably not-unique. In 
order to accomodate the possible transients in one of the above 
classes, a key factor is variability. The most variable candidates, 
BL Lacs, can vary up to a factor of ~ 100 (e.g. Mkn 501; Pian 
et al. 11998b . enough to explain the non-detections. This makes 
BL Lacs appealing candidates for possible transients. The vari- 
ability of type 2 quasars, instead, seems not sufficient. Other 
known transient source classes usually can not produce such 
high fx/fopt values. An alternative explanation could be that of 
a new kind of transient, capable of great flux variability but 
reaching much lower peak maximum luminosity. According 
to this idea, the fx/f pt would be greatly boosted by the non- 
simultaneity of X-ray and optical observations. 

6. Follow-up projects 

We wish to stress that the interpretation of BFS requires future 
multi-wavelength observations. Deeper optical and near-IR ob- 
servations should allow us to find the optical counterpart (ex- 
cept possibly for INS), and future X-ray observations should 
give spectral and temporal insights. We already started a pro- 
gram of optical follow-up. In the future, more will be known on 
these mysterious and fascinating lacunae in the X-ray/optical 
sky. 
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